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1. Introduction 

Faraday considered lines of force to be a physical substance, the basic dynamical 
object of electromagnetism, and indeed other physical theories. In his view 1 the 
particles acted on by forces were not separate entities but actually configurations of 
forces, and the lines themselves could be physically shaken, disturbances propagating 
along them with finite speed thus accounting for radiation without the need of an 
ether. 

As the quantitative description of electromagnetism was developed this point of 
view lost ground to Maxwell's which took the electric and magnetic fields, as well 
as moving charges, as the dynamical degrees of freedom 0. Gauss' law 2 in integral 
form ^> s E ■ dS = J pdV/eQ seems to support Faraday's interpretation. It measures 
the number of lines of electric force cutting a closed surface, E, as though these were 
physical objects capable of being counted and sets it equal to the enclosed charge. In 
differential form, V-E = p/e , this becomes just a differential equation that relates E 
and p, on the same footing as the other Maxwell equations that together determine 
the dynamics of the theory. From Maxwell's point of view the lines of force are no 
more than a geometric representation of the field, and the integral version of Gauss' 
law simply a statement about the ends of the curves used in that representation. 

1 "You are aware of the speculation which I sometime since uttered respecting the view of the 
nature of matter which considers its ultimate atoms as centres of force, and not as so many little 
bodies surrounded by forces.... The view which I am so bold as to put forth considers, therefore, 
radiation as a high species of vibration in the lines of force.... It endeavours to dismiss the ether, 
but not the vibrations." jjj 

2 We will use S.I. units to express Maxwell's equations H 
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The canonical quantisation of electrodynamics moves even further from Fara- 
day's picture. Based as it is on a Hamiltonian formulation, the dynamical degrees 
of freedom are the gauge potentials A , A, modulo gauge tranformations, rather 
than E and B. Their existence follows from the Maxwell equations V • B = and 
V x E = — B, which are therefore implemented as identities in the quantum theory. 
The remaining Maxwell equations, i.e. Gauss' law and 

V x B = /i J + /i e E (1.1) 

are the Euler-Lagrange equations for the theory. To quantise one has to pick a gauge. 
With the choice Aq = 0, which is most convenient for the Hamiltonian formalism, 
(|1.1|) is Hamilton's equation of motion. However the choice of gauge removes Aq as 
a dynamical variable so that Gauss' law is not recovered in this way, and has to be 
imposed as a constraint. In the quantum theory this becomes a restriction on the 
physical states of the theory. The interpretation of V • E as the generator of the 
remaining gauge transformations that preserve the gauge condition implies that the 
physical states are selected to be those that are invariant under time-independent 
gauge transformations. So Gauss' law, which from Faraday's point of view counts 
the dynamical objects in the theory, is not even valid for all the states needed to 
construct the quantum theory, but only for the physical subspace. 

There is also a conceptual difficulty with treating lines of force as dynamical 
objects. The lines, being tangent to the field, encode its direction. Its magnitude is 
represented by their density. As fields vary continuously with position they cannot 
be modeled by a whole number of lines, making obscure the concept of individual 
lines as dynamical objects. 

The purpose of this paper is to attempt to resurrect, albeit in modified form, 
the notion that lines of force can indeed be treated as dynamical objects in their 
own right. In classical electromagnetism the lines of force are fixed by the charge 
distribution, but we will consider the consequences of allowing their positions to 
fluctuate. We will overcome the conceptual difficulty by assuming that the large 
number of lines of force stretching between macroscopic charges should be treated 
using statistical mechanics, so that the classical electromagnetic field emerges as an 
average using an appropriate Boltzmann weight, and thus can vary continuously with 
position. String theory provides us with a natural identification of this weight and 
the technology to compute the average. 

In 1955 Dirac proposed a similar solution to this conceptual difficulty ||. He was 
interested in constructing a version of QED in which the electrons were created by 
gauge invariant operators which create a part of the electromagnetic field along with 
the charge. By taking this to consist of a single Faraday line of force associated with 
a quantum of flux he provided another explanation of electric charge quantisation 
(in addition to his celebrated argument based on the motion of an electric charge in 
the presence of a magnetic monopole.) In this theory closed lines of force describe 
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photons, open lines describe electron positron pairs and pair creation is described 
by the breaking of lines of force. Although the flux is supported on the line of force 
quantum superposition allows for states with continuous fields, such as the spherically 
symmetric field of a single charge. In our approach we will construct the classical 
electromagnetic field as a thermal average over a macroscopic number of lines of 
force connecting electric charges that are large in comparison with the charge of the 
electron. Dirac's theory could be taken as the microscopic description underlying 
this. 

We will begin by studying two simplified cases namely electrosatics and magne- 
tostatics. The first is easy to formulate mathematically, whilst the second already 
requires the formalism of string theory and is a useful stepping stone to constructing 
the full time- dependent electromagnetic field. By magnetostatics we mean the time 
independent magnetic field generated by constant currents flowing around closed cir- 
cuits. The lines of force associated with B form closed curves because V • B = 0, 
and Farady thought of these as dynamical objects too. However, for this particular 
example we will not focus on these lines of force. Instead we will take the dynamical 
objects that describe magnetism to be surfaces spanning the circuits the currents 
flow round. Again there is a natural weight to average over these surfaces, but we 
might expect to be impeded by the well-known difficulties encountered in trying to 
formulate sums over random surfaces that make it difficult to construct string theory 
away from its critical dimension ||. Remarkably these difficulties are absent from 
our problem even though it amounts to an off-shell calculation in non-critical string 
theory and we are able to evaluate the sum and show that it yields the Biot-Savart 
law. 

The magnetostatic problem generalises to higher dimensions in which context 
the fluctuating surfaces can be re-interpreted as the world-sheets of lines of force. 
By choosing the target-space appropriately we will show that the retarded solution 
to Maxwell's equations arises naturally as a thermal average over these lines, so 
that by invoking statistical mechanics this approach violates the usual time-reversal 
invariance of classical electromagnetism. 



2. Electrostatics 

We begin with the special case of a static electric field. Consider the field, E(x), due 
to two equal and opposite charges, ±g, placed at a and b respectively. 

. . (/ x — a ox — b . 

Ex =T^-u US-T-TI KM8- 2J 

Aireo ||x — a|| d 47reo ||x— d\\ 6 
This is the unique solution to V x E = and the differential form of Gauss' law 

V-E(x) = ^-o" 3 (x-a) - ^-o" 3 (x-b), (2.2) 
eo e 
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that decays at infinity. Briefly setting aside the former of these two differential 
equations, it is easy to see that (|2.2| ) is solved by taking the electric field to be given 
by 

E'(x) = ^ [ 5 3 (x-y)dy (2.3) 



e o Jc 

for any curve C from a to b, since for any differentiate test-function u(x) that 
vanishes at infinity: 

J V-E'(x) u(x) dV = - J E'(x)-V«(x) dV = - J 5 3 (x - y) dy^j ■ V«(x) dV 



1 [ V«(y) ■ rfy = - w(a) - ^- u(b) = / - (5 3 (x - a) - <5 3 (x - b)) u(x) dV . 

(2.4) 

( p.3|) has the same mathematical form as the Dirac string used to represent the mag- 
netic field of a monopole. In that context the position, C, of the string has no physical 
significance, and can be changed by a particular kind of gauge transformation. In 
our work, as in [Q, however, we will attach physical meaning to C, treating it as the 
position of a physical object. Now the electric field of ( |2.3p is supported on C and so 
is completely different in character to that of (J2.1|) which is supported everywhere. 



We could describe the physical object that has position C as a string of electric flux. 
To obtain ( |2.1|) we will make the assumption that the theory is stochastic in the 
sense that the positions of the flux strings are to be averaged over with a Boltzmann 
weight, e -/W , so that E(x) = (E'(x))c. The average of any functional of C, Q is 



given by the functional integral 

(n> = | f v y n 



e ~pH\y] 



(with Z a normalisation constant, so that (l)c = 1). The physical interpretation 
is that the macroscopic charge q generating the classical electric field is composed 
of many microscopic or elementary charges, of magnitude qo, each of which is the 
terminus of a line of force as in Dirac's theory. These lines are physical objects in 
thermodynamic equilibrium at temperature and there are q/qo of them. Each 
has an electric field given by ( |2.3| ) with q replaced by q , and these contributions add 
up to give the total electric field. 

We keep the end-points of all the curves to be averaged over fixed at a and b so 



the calculation (|2.4j) goes through as before for the averaged field and Gauss' law is 



satisfied. We now have to find fiH so that V x E = 0, or equivalently, so that 
x — a x — b 1 



i J Vy jj 3 ^-y)dye-^. (2.5) 



There is a natural choice for /31-L that occurs in the path-integral representation of 
quantum mechanics M and the heat-kernel connected with diffusion and Brownian 
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motion. If y(t), < t < T is a parametrisation of a path from a to b then 

,-||a-b|| 2 /(2T) 



(b|e- 



-TH \ 



a) = / Dye-/o^ 2 / 2 



(2.6) 



(2ttT) 3 / 2 

where H = p 2 /2 so 2(h\H = V 2 (b| and the eigenstates of position are normalised 
to (b|a) = <5 3 (a — b). If we take (3 = l/T and W, = Jq du (dy/du) 2 /2 then by a 
change of variable, t = Tu, we get /3T-L = J Q T dt y 2 /2. We have to set a value to T. 
This is a dimensionful quantity, and no such parameter appears in (|2.5| ) so we will 
take the limit in which T is large (in comparison to the squares of the other lengths 
in our problem.) The expectation value of the delta-function can be generated by 
functionally differentiating with respect to a source term added to /3"H: 

8 



Vy / 5 3 (x-y)c2ye-Jo d ^l 2 



c 



SA{x) 



Vye 



■fo 1 dty 2 /2+S»A( y ydy 



A=0 



(2.7) 

The functional integral inside the braces is the generalisation of ( |2.6| ) to a particle 
moving in an electro-magnetic field with vector potential iA, so the Hamiltonian Hq 
is modified to H = (p + iA) 2 /2. Thus (|2]7|) can be written as 



■(b|e- T "|a> 



<U(x 



,4=0 



Tdt ^ i '- T)H "^ ) 



-tH 



A=0 



a) 



Now when A = 0, 2 5H /(L4(x) = ip 5 3 (q — x) + 5 3 (q — x) ip, so using the resolution 
of the identity j |c) d 3 c (c| = 1, gives 



5H 



6A(x) 
so that 



.4=0 



(V c |c)) d 3 c(c| 5 3 (c-x)+ / 5 3 (c-x) |c) d 3 c V c (c| = |x) V (x| 



{21 



becomes (after setting A = 0) 

1 r T ++ 

— / ^(b|e (t - T) ^°|x) V (x|e-* Ho |a). 

2 ./n 



The normalisation constant, Z is just the right-hand-side of (|2.6|) so now we have 



1 

Z 



Vy / 5 3 (x-y)c?ye- 



(27TT) 



3/2 



C 



-b||2 



-||x-b||^ 

e 2(T-t) 



— I |a — x| 

<-> e 2i 
V 



o (27r(T-t)) 3 / 2 v (2vrt) 3 / 2 



2e 2T" 

For large T the integrand is negligible except when t ~ and t ~ T, so that in the 
limit of infinite T the integral separates into two contributions: 



dtV- 



-||a-x[[V(2t) 

2(2vrt) 3 / 2 



dt V- 



|x-b|| 2 /(2t) 



V 



1 



47r 1 1 x — a| 



+ 



1 



An\ |x — b| 



2(27Tt) 3 / 2 

which yields the right-hand-side of (j2.5|). 

When we add up the contributions to the electric field of the individual lines of 
force we obtain 



Q_ 

e o Jc 



5 3 (x-y) dy) 



Q 



x-b 



c 



47reo llx — all 3 47reo llx — bll 3 



Efx). 



(2.9) 
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3. Magnetostatics 



We now consider another simplified case, namely the magnetic field generated by a 
current that is constant in time. Although this is a different physical problem to 
that of the preceding section the mathematical description will provide a stepping- 
stone to introducing time-evolution into the description of the electromagnetic field 
generated by point charges. A constant current, /, flowing around a circuit C has 
density 

J(x)=/£5 3 (x-y)dy. (3.1) 

V • J = follows from a similar argument to ( |2.4|) . The constant magnetic field it 
generates is given by the Biot-Savart law 

B(x) = ^/ ^ X(x - y) , (3.2) 
4vr J c ||x-y|| d 

which is the unique solution to the Maxwell equations V • B = and V x B = /x J 
that vanishes at infinity. The first of these implies that the lines of magnetic force are 
closed, however we will not focus on the flux-lines, but consider the representation of 
B in terms of surfaces normal to it with density proportional to its magnitude, i.e. the 
equipotential surfaces for the (multi-valued) magnetic scalar potential proportional 
to the solid angle subtended by C at the point x. In integral form the second of this 
pair of Maxwell equations is Ampere's law 

I Bdx = / //oJ-dS. (3.3) 

JC </£' 

If C is chosen to loop around C then the right-hand-side of this is [il and the left- 
hand-side counts the number of surfaces representing B that are cut by C . These 
open surfaces all have boundary C. This is a generalisation of Gauss' law, and we 
will interpret it similarly as an indication that the theory can be rewritten in terms 
of dynamical objects which are the surfaces spanning C. In classical electromag- 
netism these surfaces are of course fixed once C is specified, but we will investigate 
the consequences of allowing the surfaces to fluctuate so that the magnetic field is 
obtained by averaging over them with an appropriate weight. 
V x B = /i J is solved by taking the magnetic field to be 

B'(x) = /i /^ ( 5 3 (x- y ) C /S(y), (3.4) 

where S is any surface spanning the current circuit C. This is readily shown by 
integrating B' against the curl of a vector test-function. Furthermore this remains 
true on averaging over all such surfaces spanning C with any weight. We will now 
endeavour to find a weight so that the averaged field also satisfies V • B = 0. As in 
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the previous example of electrostatics there is a natural choice, but before we invoke 
it we return to that used for averaging over the lines of electric force and write it in 
a form that will motivate the generalisation to surfaces. In the previous section we 
used 

(Q)c= lim \- f Vyne-tf dt * 2/2 . (3.5) 

T^oo Z J 

Consider replacing the exponent J Q T dt y 2 /2 by [|7] 

where y(£), < £ < 1 is a different parametrisation of the path and g(£) > is 
a new variable. (|3.6|) is invariant under diffeomorphisms £ — >■ £ that preserve the 
parameter interval provided that g(£) transforms as an intrinsic metric g(£)<i£ 2 — > 
#(£) ^£ 2 — #(£) ^£ 2 an d y(£) — > y(£) = y(£)- To construct this new weight we have 
to choose some value for #(£), but which value we choose will not affect the result 
for the electric field, as we will see. If we change parameter from £ to t given by 
t = T f{; y/g(?) d£' then becomes dty 2 /2 with V = T JjJ v ^c/£. Using this 
in ( p.5| ) gives the same results as T and T' tend to infinity. 

This form of the weight has a natural generalisation to the sum over surfaces 
used in Polyakov's approach to the bosonic string, 0. Let a surface E spanning C be 
parametrised by Y(£ 1 ,£ 2 ) with the 'world-sheet co-ordinates' £ a lying in some fixed 
domain D, then 

(n h = i/OTn exp (-jJL. jf>g . g ^ A ) , (3.7) 

where ^ a fe plays the role of an intrinsic metric on E, g ab is its inverse, g = det(g a b), 
and a' is a dimensionful constant. We will refrain from integrating over g a b, but 
rather choose a value for it and find that, as before, the averaging does not depend 
on the value we pick. We will now show that with this weight 

B(x) = (fiolJ 5 3 (x - Y) dS(Y) ) E (3.8) 

satisfies the Biot-Savart law ( |3.2| ). In doing so we will encounter the usual problem 
of trying to formulate string theory away from its critical dimension. 

We evaluate (|3.8|) in the standard way by first exponentiating the Y depen- 
dence, using a Fourier decomposition of the delta-function and generating dS = 
ie ab (<9Y/<9£ a ) x (<9Y/<9£ fe ) by differentiation with respect to sources: 

B(x) = ( / 5 3 (x-Y)rfS(Y)) s = / -J-^d^e ab ^x4r- [ VY e' 3 ' , 
KJ Ve J 32ttV ^ d'f dj b Z J j=o' 
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2na'S' 



,ab 



IdY dY 



zk-(x-Y) + Y-r4-|> 5 3 (^-0 ) ^ 



The dependence on k and j is separated out by writing Y as the sum of a classical 
solution to the Euler-Lagrange equations for S', Y c , and a quantum fluctuation, Y: 



Y c + Y. 



% (^99 



ab 



where on the boundary of D, Y c coincides with the current circuit C, i.e. y, and Y 
vanish. Y c can be found using the Dirichlet Green function for the Laplacian on D, 

G, 

Gd o f * -r^) 5 3 (e' - o rf? + y c (l) • 



Y c (e) 



y c (0 



Z cancels against the source-independent parts of the functional integral (including 
the functional determinants) giving 



Bfx) 



d 3 k 



^d 2 ie ah — x — e 

327r 4 a' df d} b 











2na'S" = -- <; ik + j a 4- 



ik + f 



o 



zk + j ? 



ae 6 

y c (0 + «k ■ x 



j=o 



This involves the Green function and its derivatives at co-incident points. To make 
this well-defined we introduce a regulator, e > 0, via a spectral decomposition. If 
wa(0 is an eigenfunction of the Laplacian belonging to eigenvalue A, vanishing on 
dD then we can choose it to be real and using the fact that A > take the regulated 
Green function to be H 



G(£,f) = 5>A(fl«A(f) 



A 



Let ip denote the value of this at coincident points, ip is greater than or equal to 
zero, vanishing only on the boundary dD. Furthermore, because G is symmetric we 
have that 



d_ 



G(t,i) 



l_d_ 

2d^ 
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enabling us to write S" as 



{' k+jr ^[} ' yc + ' k ' X 



It is now straightforward to integrate over k in B(x) and differentiate with respect 
to j to get 

This splits into two integrals. In the first we change variables from (£ , £ 2 ) to 77 = 
47ra'ip), supposing that on the boundary, dD, £ 2 is constant. The form of ijj may be 
found by relating it to the heat-kernel 

A 

and then using the modification of the Seeley-de Witt expansion proposed in M: 



A.TXT \ It 

r \ 

where the sum runs over all geodesic paths linking £ and £, including reflections at 
the boundary, and cr(£, £) is twice the square of the path-length. For £ = £ and 
small r the path of zero length dominates for points away from the boundary, and 
Q ~ 1 + 0(t). For points close to the boundary the shortest reflected path is also 
important, and the boundary conditions require that for this Q ~ — 1 + 0(r). So, if 
a/2 is the square of the closest distance to the boundary then 



f-(l-e-* £ >) 
Ji 4vrt V ; 



giving 



whereas for a » e 



ip ~ <r/(87re), for a « e 



ip ~ log(cr/e)/(47r) . 

Consequently, as £ moves away from the boundary rj varies from to a large positive 
value over a distance of order y/e. In the interior of the domain this large value 
suppresses the integrand, consequently we only need to consider contributions to the 
integral from points close to the boundary. As the regulator is removed, i.e. e — > 0, 
we can ignore the variation of y c with 7], replacing it with its boundary value, y, and 
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also take the 77- integration limits to be and oo. So as the cut-off is removed this 
integral becomes 



dv f de % >< (y - *) e ~ (y - x)2/ " = )r / Vl*^ > (3-9) 



J rfl 2 J 2vr 3 / 2 w ' 4vr J c ||x 

which is the Biot-Savart law. We now argue that the second integral 

d 2 £, ^ab fyc w dy_c -(y c -x) 2 /(47ra'V) 



D 



2(47r 2 a'^) 3/2 d£ a d£ 



vanishes as the cut-off is removed. Again we ignore contributions to the integral from 
the interior of the domain restricting our attention to a strip bordering the boundary 
of D and change variables to so that (|3.10|) becomes 

J r] 6 ' 1 at, 1 oa Or) 
Replacing y c by its boundary value as before gives 

where rj ranges from to h over the width of the strip. Using ip ~ a j (87re) gives 



dr) da e _ (y _ x?/v ^ 2e f h (h)_^ y _^ h ^ ^ 



This last integral is positive and less than 



2e f°° dr)_ e _ {y _ x) 2 /r] = Ae^ 



a' j rfl 2 a'|y 

which vanishes as the cut-off is removed. Consequently only the first integral (J3T 
survives and we have established that 03.81) satisfies the Biot-Savart law. 



?s/„ vuc.vr, _ 1 I d y x ( x -y) _vt^ 1 I rf y 



<F(x-Y)dS(Y)) B = — d J v || ; =Vx- — ^Z— - (3.14) 
- 47r/ c ||x-y|| 3 47r/ c ||x-y|| 

The result is independent of g a b so it is unchanged if we integrate over this metric 
degree of freedom as in Polyakov's approach to string theory ||. This independence 
is remarkable because although (|3.4j) itself does not contain g ab the computation of 
its expectation value required a regulator and the use of ip which introduce such a 
dependence. String theory calculations are replete with quantities acquiring such 
'anomalous' dependence on the world-sheet metric, which make it difficult to aver- 
age over surfaces without imposing extra conditions, for example on the dimension 
of space or the mass spectrum of excitations. When the delta-function of (513) is 
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represented as an integral over k we are effectively summing over all the spectrum, 
and yet no mass-shell condition had to be imposed restricting this sum. 

The choice of weight (|3.7p is a natural one to make in the context of string 
theory It is also natural from the point of view of magnetostatics 3 since the energy 
in the magnetic field (|3.4j) , -^^y / d 3 x B 2 , is proportional to the area of E albeit with 
a divergent coefficient and, as is well-known, the exponent in ( p.7|) reduces to the 
surface area on eliminating g a b through its Euler-Lagrange equation. 



4. Time-dependence 

We now turn to our main problem which is the representation of the electromagnetic 
field of moving charges in terms of fluctuating lines of electric flux, so that the classical 
field results from a thermal average over a distribution of these lines. (We will treat 
the lines of force by assuming that they connect equal and opposite electrical charges 
and so have finite extent. Using this approach we could treat semi-infinite lines of 
force associated with single charges by taking the limiting case in which one charge 
in each pair is sent to infinity.) If the charges of section 2 are now allowed to move 
then the electromagnetic field generated by the current-density 

/OO /' X 

8\x - a) a" dt - q 5 A (x-b)b fl dt (4.1) 
-oo J —oo 

will vary in time. Maxwell's formulation of electromagnetism is time-reversal invari- 
ant. The classical problem of computing this field by solving Maxwell's equations 
leads to solutions in terms of advanced or retarded potentials (or indeed linear combi- 
nations of the two) that have to be distinguished from each other by an application of 
common sense rather than from fundamental theory. We want to show that if instead 
the field is generated by an ensemble of lines of force in thermodynamic equilibrium 
the retarded solution arises naturally. 

Assembling E and B into the antisymmetric tensor F^ u in the usual way by 
taking F 0i = €q and /xq F^ = — allows the two Maxwell equations containing 
sources to be written as 

d*F IH , = J v , (4.2) 
which can be solved by taking F^ u to be 

-q [ 5\x-Y)dY,^(Y). (4.3) 

for any surface E that spans the world-lines of the two charges. This solution de- 
scribes a field supported on the surface E which can therefore be interpreted as the 



3 I am grateful to Benjamin Doyon for pointing this out. 



world-sheet of a line of electric flux. To obtain a field that satisfies the remaining 
Maxwell equations 

d li F vp + d v F pli + d p F^ = 0. (4.4) 
we look for a suitable Boltzmann weight with which to construct an average over S: 

F,u(x) = ~(qj - Y) d^u(Y)h ■ (4-5) 

As in section two we interpret this as an average over a distribution of many ele- 
mentary strings, each associated with basic charge go an d in thermal equilibrium, q 
rather than go enters ( |4.5|) because the individual contributions of the flux-lines must 
be summed. We will ultimately construct this average but first we digress briefly by 
considering the problem in four-dimensional Euclidean space as this provides a useful 
step towards constructing the full Minkowski space theory. Because of the similarity 
between ( |3.1|) and (|4.1|) this problem is solved by a straightforward generalisation of 
the computation of the previous section. This can be done in any dimension, but if 
we specialise to four dimensions then by taking the generalisation of ( |3.7[ ) to be 

(n h = I JvYH exp (- JL. j^ Gm • , (4.6) 



with G^v = 5^ we obtain as the generalisation of ( |3.14| ) 



da^ f db,. 



\x — a\\ 2 J \\x — b\\ 



(4.7) 



in which — a|| 2 isa Euclidean Green function for the Laplacian. 

This result can be Wick rotated to Minkowski space by x 4 — > — ix° and using an 
e prescription to encode the position of the poles. The effect is to replace the Green 
function 1 

^ 7Z ^ 2 TZn .n^ 2 I T. ( 4 ' 8 ) 



\x — a| |^ (x — a) 2 — (x° — a ) 2 + 

This does give a solution to Maxwell's equations in Minkowski space, however because 
it was obtained by Wick rotating the functional integral it has the physical inter- 
pretation of being a quantum expectation value (and therefore relevant to Dirac's 
microscopic theory jp) and not a thermal average. As a consequence the causal prop- 
erties of this solution are not those we seek, but rather they are those inherited from 
the Feynman propagator ( |4.8| ). To construct the average we want we consider the 
general problem of constructing thermal Green functions in quantum theory. (It is 
not sufficient to consider classical statistical mechanics because strings are extended 
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objects described by two-dimensional field theory on the world-sheet consequently 
this computation would be afflicted by the ultra-violet catastrophe unless we in- 
voke quantum mechanics). In the quantum theory of a dynamical variable (p with 
Hamiltonian operator H, associated eigenkets \E), and a set of time-dependent op- 
erators &i(ti), .., Cl n (t n ) the finite-temperature Green functions are thermal averages 
of time-ordered products: 

(Q n (t n ).A(ti)) T = N^e-t E (E\n n (t n )..Cl 1 (t 1 )\E) (4.9) 

E 

with N = 1/ J2e e ~ l3E an d t n > .. > t\. Making the time dependence of the operators 
explicit 

tt^tj) = e i(t ^°^%(t )e~ l ' fe ~ to)J * (4.10) 
and expressing the sum over energy eigenstates as a trace puts this into the form 

iVTr ^e-^e i( *" +1 -' o) ^e- i( '" +1 -'™ ) ^n n (t )e" i( '""*"- l) ^..fii(to)e" i( ' 1 ~' o) ^) (4.11) 

where we have introduced a smallest time to an d a greatest time t n . The functional 
integral representation can be constructed in the usual way to give 

i j V<pe i So*>*n n (t n )M 1 {t 1 ) (4.12) 

where L is the Lagrangian related to H . The contour C consists of three straight- 
line segements. The first, C\ runs just below the real axis from to to t n+ \ — ie' . The 
second, C\ runs from t n+ i — ie' to to — 2ie' and the final segment, C3 runs vertically 
down from £ — 2ze' to to ~~ 0- e ' is included to ensure convergence in the exponential 
factors. Note that the Q are inserted only on the first segment, G\. The trace is 
computed by identifying the values of the integration variable tp at to and to — ifi. 
(See pi and references therein for a discussion of complex time contours and thermal 
Green functions). 

We will interpret this construction as specifying the target space for our finite 
temperature theory. The world-sheets of the electric flux lines should wrap around 
the contour C. We will take /3 to be small (on the scale of the typical distances 
|x — a| and |x — b| involved) so that we can neglect the contribution to that part 
of the world-sheet on C 3 . Ultimately we should send x — >■ —00 and x n+ i — > 00 
resulting in two infinite segements, C\ and C2. The edges of the sheet are the world- 
lines of the two charges, and these are duplicated on the two segments, however 
the operator whose Green function we are computing is restricted to C%. There are 
thus two contributions to the thermal average. The first comes from the part of the 
world-sheet on C% and gives the Minkowski space version of ( |4.7| ) based on the Green 
function (|4.8|). The second comes from the part of the world-sheet on Ci- Because 
x° G C\ but a e C2 the Green function is 

(x - a) 2 - (x° - a° + ie") 2 
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instead of ( |4.8|) . (e" is positive, proportional to e' and has an irrelevant dependence 
on x° and a ). When the direction in which C2 is traversed is taken into account these 
two contributions combine to give an expression like ( f4.7| ) but with the replacement 



\x — a\ 



12 

1 1 



(x - a) 2 - (x° - a ) 2 + ie (x - a) 2 - (x° - a ) 2 - 2t{x° - a°)e" ' 

(4.14) 

Now this is a representation of the retarded Green function 

%iS (x° — a — |x — a|) 



-2me(x° - a ) 5 ({x° - a ) 2 - (x - a) 2 ) 



x — a 



Taking into account the factor of i that 5 4 (x — Y) acquires under Wick rotation we 
find that the contribution to the Minkowski space thermal average 

-( / 5\x-Y)d^, u (Y))^ (4.15) 

coming from the edges of the world-sheet that are the world-lines of the charges is 
n / f , 5(x°-a - |x-a|) f „ 5 (x° - b° - |x - b|) 

d v / da » rn i / dh » rn u 

V/ 47r x — a / 47r x — b 



(4.16) 



, . , 8(x° -a -|x-a|) f „ 5 (x° - b° - |x - b|) 

— dy I / da u — : — — / db u — : — - 

47tx — a / 47tx — b 



Since we are taking small so as to be able to ignore C3 the trace is evaluated by 
identifying the string configurations at the end of C\ and C2 at times t and to — 2ie' 
respectively, and then integrating over them. Before the integration is done the 
configuration of the string constitutes two fixed edges of the world-sheet described 
by two copies of the same curve displaced through 2ie' in time. Taking the edge 
attached to the end of C\ to be = c M , say, there is a contribution to (|4.15|) of 



4, ( , de /^ :l-'r c|) N ) -m -» i . (4.17) 



47r|x — c| / \J 47r|x 

Because this is supported on null-rays through Y^ = c M it vanishes for fixed x when 
we send to — > — 00 (provided that the world-lines of the two charges, and 6^ 
are not themselves null). Consequently only (|4.16[) survives in this limit. Summing 



the contributions from all the elementary charges go by simply multiplying this by 
q results in the retarded solution to the Minkowski space Maxwell equations, F, 



ret 
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so we arrive at the conclusion that the retarded solution arises naturally from the 
thermal average: 



in the limit of small /3. Notice that as in section two we have taken the temperature 
to be large. 

5. Conclusions 

The search for a unified description of physical phenomena is an old one. Faraday's 
idea that the basic atoms of such a description are lines of force may appear as 
no more than a quaint notion that was rapidly side-lined by Maxwell's quantitative 
description of electromagnetism in terms of fields. However, from a modern perspec- 
tive informed by string theory, his notion seems remarkably prescient, as indeed do 
Dirac's related ideas of over half a century ago. 

We have used some of the world-sheet technology of string theory to show that 
classical electromagnetism can be interpreted as a consequence of the statistical me- 
chanics of lines of force. This provides an underlying stochastic description of what 
is usually taken to be a deterministic problem encoded in the partial differential 
equations of Maxwell's theory. The result of this is that the retarded solution arises 
naturally, breaking the time-reversal invariance of the classical theory as a conse- 
quence of invoking thermodynamics. The direction of time is picked out because 
when we invoke the statistical mechanics of flux-lines we necessarily have to use 
quantum theory to avoid the ultra-violet catastrophe associated with the statisti- 
cal mechanics of fields, since the world-sheets of strings are effectively field theories 
with one space and one time dimension. Quantum theory requires a Hamiltonian 
operator with a spectrum that is bounded from below. This allowed us to introduce 
e insertions as convergence factors in our discussion of time dependence in section 
four, and it is these that pick out the retarded, rather than the advanced, solution 
to Maxwell's equations. 

Electromagnetism, and gauge theories generally, are intimately connected with 
string theory and arise in a number of ways so we should consider how our compu- 
tations are related to these standard approaches. 

Conventionally photon vertex operators are open-string operators inserted on the 
world-sheet boundary and correspond to coupling to an external field by adding 
to the action § dy^ so that A^ acts as a source for the operator. The scattering 
amplitudes of photons of definite momenta and polarisation are thus obtained 
from the expectation values of vertex operators q § e tk ' y E ■ dy on the world-sheet 
boundary, and the mass-shell condition k 2 = results from the requirement of the 
decoupling of the scale of the world-sheet metric. In our flux string picture the 
electromagnetic field strength is represented by a closed-string vertex operator 




(4.18) 
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that is inserted on the interior of the world-sheet and corresponds to coupling to an 
external field by adding J qB^ v dE^ to the action. We can make a connection 
between the conventional open string vertex operators and our closed string ones by 
recalling the LSZ formalism in QED. This gives scattering amplitudes as time-ordered 
expectation values of on-shell field operators 

J d 4 xe tk - x E"d 2 A^x) = - J d^xe^E^F^x) (5.1) 

since in QED F^ v = d^A v — d u A^ is an identity (although it only holds on average 
in our work), and k ■ E = 0. If we were to calculate this in terms of the flux string 
picture using (|4.5|) we would calculate the expectation value of 

q [ d 4 xe tk - x E"d' / [ 5\x - Y) dE^(Y) = q (f e lk ' v E ■ dy (5.2) 

J JT, JdY, 

which is the usual open string insertion. 

The averages over the flux string configurations that we have needed to construct 
F^y can be computed in any dimension because the world-sheet metric decouples from 
the calculation, consequently there is no mass-shell condition and we are free to use 
this off-shell 5-function insertion. So our calculation is essentially one in non-critical 
string theory. 

Gauge theories also arise from strings stretching between coincident D-branes. 
Our flux-lines stretch between the world-lines of electric charges which may be con- 
sidered as DO-branes, however they are not coincident so that the flux-lines would 
not be associated with massless excitations. 

There is a further way in which our flux strings differ from usual string theory, 
and that is in their interactions. We have computed the classical electromagnetic field 
from the statistical mechanics of large numbers of flux-lines, but we could instead use 
this approach to pursue further Dirac's idea of building Quantum Electrodynamics 
from elementary flux lines. In the Euclidean functional integral approach to the 
quantum theory of the electromagnetic field coupled to charged particles we need to 
compute the expectation value of operators by integrating over the gauge field 

(Q 1 ..Q n ) A = ^-jvAe~y F ^ F ^ d4x Q 1 (A)-^n(A). (5.3) 

It is sufficient to take the operators to be Wilson loops, Jill, i- e - 

ilj = exp — iqj (p A^ dx 11 , 

for arbitrary Cj. We will not address the dynamics of the charges, which requires 
integrating over the Cj, but focus on eliminating A^. The integral over the gauge- 
potential can be computed exactly as a functional of the curves Cj, because the 
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exponent is quadratic in A^, and the result is 




where F° v is the classical (Euclidean) electromagnetic field generated by the charge 
density iqj § c 5 4 (x — y) dy^ , and Sj is any surface spanning Cj. We now represent 
the classical electromagnetic field as an average over flux-strings to obtain 

(n 1 ..n B >A = ex P x;^( I dnx)%5\x -Y)dnY)Th^ (5.4) 

i,3 

where we have introduced an extra averaging over the Ej to obtain a more symmetri- 
cal result. Although this is not a complete theory it is clear that the basic interaction 
between the flux strings is a contact interaction rather than the splitting and joining 
interaction of conventional open string theory. 
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